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We develop a relativistic Coulomb-corrected strong field approximation (SFA) for the investigation of spin 
effects at above-threshold ionization in relativistically strong laser fields with highly charged hydrogen-like 
ions. The Coulomb-corrected SFA is based on the relativistic eikonal-Volkov wave function describing the 
ionized electron laser-driven continuum dynamics disturbed by the Coulomb field of the ionic core. The SFA 
in different partitions of the total Hamiltonian is considered. The formalism is applied for direct ionization of 
a hydrogen-like system in a strong linearly polarized laser field. The differential and total ionization rates are 
calculated analytically. The relativistic analogue of the Perelomov-Popov-Terent'ev ionization rate is retrieved 
within the SFA technique. The physical relevance of the SFA in different partitions is discussed. 

PACS numbers: 32.80.Rm,42.65.-k 



I. INTRODUCTION 

Since the pioneering experiment by Moore et al. ||T| laser 
fields of relativistic intensities (exceeding 10'** W/cm^ at an 
infrared wavelength) have been applied for the investigation 
of strong field ionization dynamics of highly charged ions 
l^-lHlj see also ID. A lot of effort has been devoted to the 
numerical investigation of the dynamics of highly-charged ions 
in a super-strong fields |10-20|. 

A common analytical approach for strong field atomic pro- 
cesses is the strong field approximation (SFA) ll2TH23i which 
has been applied for the treatment of relativistic effects Il24l|25l . 
The main deficiency of the standard SFA is that the influence 
of the Coulomb field of the atomic core is neglected for the 
electron dynamics in the continuum and the latter is described 
by the Volkov wave function f26l. While this is well-justified 
for the ionization of a negative ion, for atoms and, moreover, 
for highly-charged ions it is not valid and the standard SFA 
can provide only qualitatively correct results. A modification 
of the theory is required to take into account the influence of 
the atomic core on the free electron motion. 

The Coulomb field effects during ionization have been 
treated in the Perelomov-Popov-Terent'ev (PPT) theory based 
on the imaginary time method 11271 12B1 . In the PPT theory the 
barrier formed by the laser and atomic field is assumed to be 
quasi-static and the tunneling through it is calculated using the 
WKB -approximation. The relativistic PPT-theory ||29l - l32l has 
provided the total tunneling rate. However, spin effects are not 
investigated thoroughly ll33l . 

The standard SFA technique has been modified to include 
the Coulomb field effects of the atomic core. A heuristic 
Coulomb- Volkov ansatz has been used for this purpose, see 
||34| -|45 1. In a more rigorous way, the SFA is modified by re- 
placing the Volkov wave-function in the SFA transition matrix 
element by the, so-called, eikonal-Volkov one 1461 i47l . The lat- 
ter describes the electron continuum dynamics in the laser and 
the Coulomb field. Here the laser field is taken into account 



exactly, while the Coulomb field via the eikonal approxima- 
tion. The nonrelativistic Coulomb-coiTected SFA based on 
the eikonal-Volkov wave function has been applied recently 
for molecular high-order harmonic generation |48 -50|. The 
eikonal approximation has been generalized to include quan- 
tum recoil effects ll5ni52l . The relativistic SFA based on the 
generalized eikonal wave function has been proposed in f53l . 
However, the final results have been obtained only in the Born 
approximation. 

In this paper we develop the Coulomb corrected SFA for the 
relativistic regime based on the Dirac equation, extending the 
corresponding nonrelativistic theory, see paper I (the first paper 
of this sequel), into the relativistic domain. This will allow us 
to calculate spin-resolved ionization probabilities taking into 
account accurately the Coulomb field effects of the ionic core. 
In the Coulomb corrected SFA, the eikonal-Volkov wave func- 
tion is employed for the description of the final state instead 
of the Volkov one. The influence of the Coulomb potential of 
the atomic core on the ionized electron continuum dynamics 
is taken into account via the eikonal approximation. Direct 
ionization of a hydrogen-like system in a strong linearly polar- 
ized laser field is considered. Two versions of the relativistic 
Coulomb corrected SFA are proposed that are based on the 
usage of different partitions of the total Hamiltonian in the 
SFA formalism. The physical relevance of the two versions is 
discussed. 

The plan of the paper is the following: In Sec.|ll]the rela- 
tivistic Coulomb-corrected SFA is developed. The differential 
and total ionization rates for hydrogen-like systems are derived. 
The modified SFA using a specific partition of the Hamiltonian 
is considered in Sec. |lll] and is used for the calculation of 
ionization rates. 
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We consider the interaction of a highly-charge ion with a 
strong laser field in the relativistic regime which is described 



2 



by the time-dependent Dirac equation: 



1 V*^"^' 
c c 



(1) 



where y^' are the Dirac matrices, V*''^ = -/c/r is the Coulomb 
potential of the atomic core, k the typical electron momen- 
tum in the bound state, determined by the bound state energy 
via - Ip - Vc^ - c^K^, Ip the ionization potential, is 
the 4-vector potential of the laser field (atomic units are used 
throughout). The transition amplitude for the laser induced 
ionization process from the initial bound state cfr. into an exact 
continuum state t// with an asymptotic momentum p can be 
written 1241 



i/ 



d'^x^ir, t)A^f(r, t) 



(2) 



with 4 = y^A^i. Equation ^ is still exact with the exact wave 
function i/'(r, f) describing the dynamics of the electron in the 
ionic and the laser field. Here, we have used the standard 
partition of the total Hamiltonian: 



H - Ho + Hint, 

Ho = ca ■ p + pc^ + V^'\r) 



(3) 
(4) 
(5) 



where a - y^y, fi - y" are the Dirac matrices, and c the speed 
of light. We consider a hydrogen-like highly charged ion in 
the relativistic parameter regime, i.e., the wave function of 
the initial bound state 0: fulfills the Dirac equation with the 
Hamiltonian //q and is given by ll54l 



where the bispinor 



^3/2 



-ilKrY 



X exp |^-/fr - - 7^)? j v,-, 



Xi 



(6) 



(7) 



describes the spin-up and -down states (/ = +), with the two- 
component spinors;if+ = (1,0) and;^^^ = (0, 1), respectively, 
and the Pauli-matrices cr. 



The hydrogen-like system interacts with a strong linearly 
polarized laser field 



£(77) - -Eo cos(w?7), 



(8) 



where?; - k^'x^jio, andA:'^ - ((y/c,k) is the laser 4-wave-vector. 
[the coordinate and the momentum projections are defined as 
r£ = r ■ e, rj; = r ■ k, = p ■ e, = p ■ k, and = P ■ (k X e), 
with the unit vectors k and e in the laser propagation and the 
polarization direction, respectively]. The vector-potential is 
chosen in the Goppert-Mayer gauge: 



A^' = {if, cA) = (-r ■ E, -k (r ■ E)). (9) 
Generally, the SFA in different gauges do not coincide and 
correspond to different physical approximations [55-57|. In 
paper I we have seen that in the nonrelativistic regime the 
Coulomb-corrected SFA in the length gauge leads, first, to 
rather simple expressions for the Coulomb-corrected ionization 
amplitude, see Eq. (1.28) [refers to Eq. (28) of paper I], which 
is due to the cancellation of the r ■ E interaction Hamiltonian 
in the matrix element by the Coulomb-correction factor; and, 
second, to results coinciding with the PPT-ionization rates. 
As the latter provides a good approximation to experimental 
observations, in the relativistic regime we choose a gauge 
which generalizes the length gauge into the relativistic domain 
ll58l . that is the Goppert-Mayer gauge defined by Eq. (|9]l. 

In the conventional SFA the exact continuum state is ap- 
proximated by the VoUcov-wavefunction which is identical to 
the first order WKB-approximation of the continuum electron 
in the laser field. A systematic improvement of this approxi- 
mation is achieved in the relativistic Coulomb-corrected SFA 
by the replacement of the exact wave function if/ir, t) with 
the relativistic eikonal-Volkov wave function. Similar to the 
non-relativistic case, see paper I, first we apply the WKB ap- 
proximation to the solution of the Dirac equation ([TJ, then the 
resulting equations will be solved via a perturbative expansion 
in the Coulomb potential V''*. 

Due to the gauge covariance of the Dirac equation, we 
switch to the velocity gauge with its vector potential A^' 



^ drj^ir] ) and after solving the 



Af'iT]) = (0, cA), A 
wave-equation go back to the Goppert-Mayer gauge, 
quadratic Dirac equation in velocity gauge becomes 



The 



ihy'^d^ + -yAf, - /— + c\ Unyi'd^ + -yfAf, - /— - c) i/r = 0, 



2a2 



iti 



p y(c)2 

f +2A-d + a- vy^') - IV^'^do) +-^ + — 



i/r = 0. 



(10) 
(11) 



where A' = ^ and we have inserted Tl to indicate the WKB-expansion. Let us assume that the solution has the form ip - e"'". 
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then the corresponding equation for S is 



A" 

C C / \ C 

n 

The WKB expansion S - So + -S i + . . . yields following equations up the first order in Ji/i 

I 

M". s,S^S,-^i^.K-^''-'-^-^. (12, 



I I ' c c 



^\ if^^:^': 2A^^5i .2.^^^' 2V^^^ a^^V^ 

- : IOuSqctSi =-oSq+ 1 . (13) 

I c c " 



These equations can be solved perturb atively with respect to 
the potential term V^'^K If we define 



So 
Si 



c(0) , r.(l) 

sf^ + s\'\ 



the following equations can be derived 



5,5 J,' V 



(0) 



id.sf^" ^d^sl'^-H', 

a2c(l) 



(14) 
(15) 



(16) 

(17) 
(18) 



(19) 



where n^' 



-d^S 



(0) 



( 16 1 gives the Volkov-action, i.e.. 



H is the relativistic momentum. Eq. 

c 



.(0) 



-p ■ X ■ 



c(k- p) 



A-p + ^]dT]'. 



(20) 



(21) 



On the other hand, the solution of Eq. ( 18 1 is 

(0) _ 

1 " 2c{k-p)' 

Here we should note that S and S j*'^ generate the full Volkov 
solution. The corresponding equations in the Goppert-Mayer 
gauge can be found via a gauge transformation with an gener- 
ating function = ^(77) ■ r. After a coordinate transformation 
from (f, r) to {t], r) they become 



5[,V77) = (p + A(77)--k|T 



E + 



J 

(p + A(7,')/2)-A(7,') 



A 



c(0)/ N (1 +k-Q')(Q'-A) 

^1 ^ 2^A ' 



(22) 



(23) 



with s - c -sjc^ + and the constant of motion A = k ■ p/a> 
Further, Eqs. ( 17 1 and ( 19 1 can be solved via the method of 
characteristics as 



d<T 



dS 



(1) 



d(T 



Vic 



,2c(l) 



+ 



2c 2 
where the trajectory is given by 

dxi" _ 
dcr 



(24) 
(25) 



(26) 



This leads to dcr - drj/A, with the relativistic kinetic momen- 
tum in the laser field 



nil]) = p(T])^VSf' +A(i]) 

,(p + A(;7)/2)-A(77) 



p + A{t]) + k- 



cA 



and the corresponding relativistic kinetic energy 



C7T^(7]) = sir]) = -d,Sf^ = e + 



(p + A(7;)/2)-A(?7) 
A 



Then the solutions of Eqs. (24 1 and (25 1 read 



f 



drf 



a ■ Vy(") (r(?7')) 



(27) 



(28) 



(29) 



V^'Hr(Ji'))doS'; 



2c 

(0) 



(30) 



|2c(l) 



■d,S'^>d^S^ 



(0) 



d^S 







with the relativistic trajectory of the electron in the laser field 
^ drj p(t] )/A, starting at the ionization phase t] 
with coordinate r. 

We can evaluate the explicit parameters for the applied 
eikonal approximation estimating the imposed conditions for 
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the expansion of Eqs. 

5f « S?' 



14b and 



si'' « 5'°' 



and 



'0 ^ "0 

J ^-^ ^ g . For this purpose we use order of magnitude 
estimations: 



(0) 


(1) 


,(0) 

-(1) 





(31) 
(32) 

(33) 
(34) 



Here, the first two equations were already derived in paper 
I, see Eq. (1.22), £„ = {2Ipfl^. Further, taking into account 
the typical time during tunneling Tc ~ y/co = k/Ei), with the 
Keldysh parameter y - cuk/Eq ETII . the value of the typical 
velocity of the electron during tunneling f-Ec ~ k [k ~ -^J^Ip], 
the interaction length rEc ~ ^EcStc, the uncertainty of the 
initial time 6tc ~ 1/ ^/kEq [from the saddle-point integration 
St^ ~ 1/5 (fj)], s{ri)c ~ c^, A ~ 1 [the electron is at rest at 
the tunnel exit], and A ~ EqTc, we come to the following 
conditions for the applicability of the eikonal approximation: 



Thus, the relativistic eikonal- Volkov wave function reads 
(A/(r,f) 

(1 + k ■ or) or ■ A(77) -I- drj'a ■ VV''^ (r{T]')) 



1 + 



2cA 

''"■^ exp[i5[,°\r,77) + /5[;'(r,7,)]. 



V(2^ 



(42) 



with the bispinor Uf for the spin-up and -down final states 

(/ = ±) 



Uf 



[ ^]{\+slc^)l2xf\ 



(43) 



where = (1,0) x- — (0, 1). The last term in the pre- 
exponential [oc VV*''^] describes the spin-orbit coupling during 
ionization (under-the-barrier motion). Via a p/c-expansion of 
the Dirac equation in the atomic and the laser field it can be 
shown that the spin-orbit coupling, given in the Hamiltonian 
by the term 



H: 



so 



■■[p(//)xVyW(r(//))l 

4c^ ' 



(44) 



« 1 



and 



^ 1. 



(35) 



The action function 5 j'* can be approximated by the first term 
in Eq. dSOli: 



Ji; 



, a ■ Vy("' (r(?7')) 
2cA ■ 



(36) 



In fact, the order of magnitude of the terms in r.h.s. of Eq. pO\ 
are 




(37) 
(38) 
(39) 



and Tn is vanishing because in the tunneling region AV*'^^ = 0. 
From the latter, we estimate the ratios: 




can be evaluated along the most probable trajectory as follows 

\3/2 



Pk <5y(^> PE 5y<^> 

n-SOT^c Tc, —:r— Te 

ove ark 



(45) 



where the typical values for pe ~ k, pk ~ k^/c, ~ rEcK/c 
have been used. It is of higher order smallness than the terms 
kept in the eikonal wave function and is neglected in the further 
calculation. 

Finally, the ionization amplitude in the Coulomb-corrected 
SFA will have the following form in the relativistic regime: 



cx 



xexp[-/5(77)]. 



(46) 



with the spatial parts of the initial |0r, i) and the final |p(7?)r, /) 
spinors and the contracted action 



5(7?) 



f 



dr( 



[p+A(?70/2]-A(?70 
A 



.(47) 



In the adiabatic regime (w <K Up, with the pondero- 
motive potential Up — E^Aco^) the time integration in the 



« 1, 



(40) 
(41) 



amplitude of Eq. (46 1 is calculated with SPM. As in the non- 
relativistic case, here also the disturbance of the saddle-point 
conditions due to the Coulomb field is neglected. This is in ac- 
cordance with our approach to take into account the Coulomb 
field perturbatively in the phase of the WKB wave function: 
therefore, only the first integrand term may be maintained '^o * is already proportional to the Coulomb field V^^Xr), there- 
fore, in the trajectory r(;7) it should be neglected. Mathe- 
matically this is justified as one can see from the estimation 



in Eq. (30i. Further, one can see from Eqs. (33i and (34i 
that 1 and 5*" < 

corresponding exponents 



that «: 1 and S'j'' <K 1, which allows us to expand the 



Ip yjEolEa [similar to estimations in Eqs. 
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It shows that is negligibly small compared to Ip, and, 

consequently, the saddle-point condition d,^S {rjs) = is not 
disturbed. The latter leads to the following two saddle points 
per cycle: 



Pe 2A{e - + Ip) - pI 



Eq/oj 



(48) 



tJT]2 - n - arcsm 



Pe 



2Aie-c^+Ip)-pl 



EqIu V (£o/w)2 



The amplitude is now evaluated at these phases. First, we 
consider the action S {rjs) in the exponent: 



Im {§(!],)} ^ 



[2A(s - + Ip) - pI] 
3|E(77o)|A 



3/2 



(49) 



with |E(77o)| - Eo^Jl - (pe/(Eo/u))^. Since the real part gives 
an unimportant phase in the resulting amplitude, only the imag- 
inary part is given. This function in the exponent dominates 
the momentum distribution and determines the maximum of 
the momentum distribution which is located at a parabola with 
Pk - Ipl'ic + /?|/2c(l -I- Ip/3c^), pb = 0. In the following, we, 
therefore, evaluate the less important pre-exponential factor 
only at this parabola and neglect deviations from it. In evalu- 
ating at the saddle point rj - rjs, note that the integration 

in 5q ^ starts at the saddle point rjs, when the electron enters 
the barrier due to the effective potential, and goes to infinity 
when the electron is far away from the core. However, the 
upper integration limit could be set at the moment 770, when 
the electron leaves the barrier since further integration over 
the continuum leads eventually only to an unimportant phase 
of the amplitude. The time which the electron spends under 
the barrier is short compared with the laser period and the 
integrand can be expanded around the saddle point 77,: 



dv'j^ (50) 
e(Tj,) + 5,e(77,)(77' - 77.) + ^^(/?' - r],f 



r + E^(7,' - 7?,) + ^(77' - 77,)2 + -^(77 - 77.)3 



where ¥i{ri) - -£(77) - k(p(77) ■ E(77))/cA is the Lorentz-force 
due to the laser field. As the correction factor is evaluated 
at the parabola with pi; - Ip/3c + /?|/2c(l + Ip/3c^) and 
Pb - 0, we derive the trajectories at these special values of 
momentum. During the motion under the barrier the magnitude 
of the coordinate variation in the k-direction (imaginary) is 
significantly smaller than in e-direction [rk{r])/rE{ri) ~ k/c] 
and the integrand can be expanded by rtiri). It is taken into 
account also that the ionization event happens close to the 
laser polarization axis: rj^ <k r^. The integral is, therefore. 



approximated by: 



-iS 



"(r,77.)= dri' 
J lis 

e(rj,) + 5,£(77,)(77' - 77,) + ^^^(77' - 77,)^ 



IK 

A? 



be + 



drj' 



A 

iK 



'-^(^'-ns)'\ 



E(n,) 

2A 



,(51) 



where 



i^pTprE 2Ip 7£(77o)y27p 

rkiri) = z — (?7 - r]s) + z (1 - rjs) 

3c ic 2c 



(52) 



Psiris) = i 



e(T]s) 

d,,sij]s) 
a2 



C^-In 



^2A(s -c^ + Ip) - pI ^ i^p[l - 5Ip/36c^) . 

(53) 
(54) 

7.£(77,)£(77o)/A (55) 
df,^^s(jj,) = £(77o)VA (56) 

Other parameters are 07770 - arcsin [-pE /E[)/cj], k « ^J2Ip( 1 - 
/p/4c2), A ^ 1 - Ip/3c\ rBiii) ^ 0, £(77,) ^ £(770). To 
have compact expressions, expansions on the small param- 
eter Iplc^ have been used above [Ipjc^ ~ 0.25 for hydrogen- 
like uranium]. The starting coordinate of the trajectory 
r is found from the saddle-point condition of the integral 
Jii''rexp[-/p(?7i) ■ r - Kr], that is rj/r, = p(77j)/(;7c). With 
these parameters the integral in Eq. (jSTJ can be calculated: 

Qr = exp[-/5<'>(r,;7,)] 



exp 



exp 



21 p 


1 - 


6c- 




c2 


A'- 


ip 




21 p 

c2 






(57) 



with the small parameter A--r- E(T]s)/4Ip. In the expression 
above only the leading order term in l/A are retained. For 
justification of our approximations, we compare in Fig.[T]our 
analytically derived Coulomb-correction factor Q,- with the 



exact result of a numerical calculation of Eq. ( 50 1, and with 
the Coulomb-correction factor of the relativistic PPT ||3T|| . De- 
viation of our approximate Coulomb-correction factor from 
the exact one occurs mostly due to a linearization of the A- 
dependence in the analytical expressions which, however, is 
needed for the further analytical integration. 



Using the Coulomb-correction factor Q,- of Eq. ( 57 1, the 



spatial integration in Eq. (46 1 can be carried out. The Coulomb- 
corrected pre-exponential matrix element reads 

mfi(ii.s) = {p(T]s)r,f\r ■ E(77,)(l - k ■ a)Qr\Or, /). (58) 



It is remarkable that the Coulomb-correction factor in Eq. ( 58 1 
cancels the dependence of the matrix element on the elec- 
tric dipole operator r ■ E and, approximately, also the pre- 
exponential term r^'p/'^' of the initial state wave function [when 
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FIG. 1: The Coulomb-correction factor vs the parameter A for 
y = 0.01 and Ip/c^ = 0.25 (equivalent to Z = 91): (black, solid) the 
numerical calculation via Eq. ( |50[ l without approximations, (blue, 
short-dashed) the analytical result of Eq. ([57j, (red, long-dashed) the 
Coulomb-correction factor of the relativistic PPT iSli . 



the approximation r ^ rE, i.e., ionization happens at the laser 
polarization axis, is used]. This is a consequence of the ap- 
plied Goppert-Mayer gauge which significantly simplifies the 
calculation of the ionization matrix element. 

In this paper we are concerned with spin-unresolved ion- 
ization probabilities. Therefore, we are free to choose the 
orientation of the quantization axis. In the following we as- 
sume that the initial spinor of the bound state is aligned along 
the laser magnetic field direction. The spin-dependent matrix 



element in Eq. (58 i in this case yields: 

r , 2c 

mfi(i],) = M r— 
J n- 



2-k 



8Kh 



2V2^ \ r(3- ^4) \\EiTls 



X exp 



2I„ 



exp [-/p(?7,) -r-Kr] 



(59) 



2'--(2y4-i?|£(;,,)pexp[|^] 



where the factor m/, for different spin transitions is 
V2c(2c + i>£)(l + f) 



mfi, (60) 



(61) 



/?2c(l68c4 - + 142cV| - 8/c/4 + 25pl) 

24 V2(2c - ipE) (2c2 + plf'^ ^4c2 -h pi 
V2c(2c-/p£)(l-f) 



(62) 



^8c4 + 6c^pI + p\ 
/?2c(l68c'* + + U2c^pI + 8icp\ + 25 p^) 

24 V2(2c + ipE) (2c2 + plf^ ^4c2 + pi 



with j3 - ^2Iplc and m_+ = m+_ = 0. Eqs. (61 1 and 



62 1 are 



As in 



exact in peIc, the PeIc term can be important at ^ » 
the non-relativistic case, see paper I, the matrix element has a 
singularity at the saddle point 77, and the modified SPM has to 
be apphed which induces an additional pre-exponential factor 



1 



1 



[p(77,)2 + /.2]2 4 (p(77,) ■ p(;7,))2 (77 - ,7,)2 



(63) 



The total pre-exponential factor after the 77-integration in 
Eq. ([46]) becomes: 



with 



V^7r(3-^)|i<(77,)p 

V2c(2c + i>£)(l +f) 
^8c4 + 6c2p|+p4 

/3^c(20c* + 2ic^pE + Uc^pI + icpl + 2pl) 
+ - 

3 V2(2c - ipE) (2c2 + plf^ ^4c2 + pI 
^/2c(2c-ipE){l-i) 
^Sc' + 6c^pI+p1 
/3^c(2Qc* - 2ic^pE + Uc^pI - icpl + '^Pe) 

H -■ 

3 V2(2c + ipE) (2c2 + plf^ ^4c2 + pI 



(64) 



(65) 



(66) 



and = m+_ = 0. In the last step the following expansion 
is used 



2mS%ri,) 



V2^ 



4 (p(77.) ■ P('7.,))' 4(2/,)3/4|£(;7,)|3/2 \ 24c2 



1 + 



41/„ 



(67) 



We note that an arbitrary configuration of the quantization 
axis can be accomplished via a rotation. Explicitly, since both 
the initial and the final state have two independent solutions 
"up" and "down", i.e., nij - +1/2 for j - 1/2, the rotated state 
is given by 



|l/2,m>' = ^|l/2,m')D„,,,„(7?) 



(68) 



where D,„' „,{R) is the j - 1/2 representation of the rotation 
matrix. Since the initial quantization axis is along the laser 
magnetic field, any quantization axis can be aligned with two 
angles, see Fig. |2] with the rotation matrix 



-'■^'/2 cos(f2/2) -e-'^'^2 sin(^2/2)^ 
'■^'/2sin(^2/2) e'^'/2cos(^2/2) . 



(69) 



Therefore the transition matrix element m'j. for arbitrarily 
rotated spinors is given by 



DjjnijkDki 



(70) 



FIG. 2: Any quantization axis can be aligned witii two angles. 



in terms of the initial one nifj. For instance, in the case the 
quantization axis is along the laser electric field the matrix is 



2 



.U.\ 

2 



1+1 1+/ 

2 2 



(71) 



while in the case the quantization axis lays along the laser 
propagation direction, it yields 



Dfi 



( I 




V2 




I 

^ V2 


1 

V2 > 



(72) 



and the calculation of the corresponding transition amplitudes 
is straightforward. 

The differential ionization rate averaged by the initial spin 
polarization and summed up over the final polarizations is 
defined as 



dw 



(l'w++('7.)P + |ot+-('7.)P + |ot-+('?.)P 
+\m^-(j]sf) exp{-2 Im [^"(t;,)]), (73) 
which with the help of Eq. (|64]) will read 

-(i-i-it-^)(2/,)^('-:^).exp(it) 



xl + 



;r2r(3-^)|£(;7..)r^(l + i)' 

-2 [2A(£-c2 + /^)-p|f' ^ 



12c2 



exp < 



3|E(77,)|A 



(74) 



Note that there is an additional factor of two, that arises due 
to summing up over the two saddle points per laser cycle. The 
momentum distribution at ionization in the relativistic regime 
is plotted in Fig. [3] As it is already mentioned, the distribution 
is located around a parabola, see also Il59]j63l . 

Characteristic features of the momentum distribution in the 
relativistic regime are the large parabolic wings corresponding 
to the large longitudinal momentum (along the laser propaga- 
tion direction) which arise during the free electron motion in 
the laser field. However, a more interesting feature of this dis- 
tribution is the small shift of the peak of the distribution with 




FIG. 3: The differential ionization rate for the parameters Ipjc' = 0.25, 
and Eg/Ea = 1/25: (a) relativistic calculation via Eq. {74\ ; (b) non- 
relativistic calculation via Eq. (1.31), Ay = ^4EJEq{EoI(jj) is the 
longitudinal momentum width. 



respect to the nonrelativistic one. While the nonrelativistic dis- 
tribution peak is at pz -Q, pk - 0, in the relativistic case it is 
shifted to pe = 0, pk- Ipl^c. As we have shown in |64|, this 
momentum shift is due to the under-the barrier dynamics and 
arises already at the ionization tunnel exit but, notwithstanding 
of the large momentum transfer from the laser field during 
the electron excursion after the ionization, the characteristic 
momentum shift of the peak of the distribution is detectable 
far away from the interaction zone at the detector. 



When the exponent in the differential ionization rate is ex- 
panded quadratically around the parabola pE - Pe, Pb - ^Pb, 
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Pk - /p/3c + /?|/2c(l +/p/3c^) + 5/7/i, the rate expression reads 
^ 2--^(l.||-i-^)exp(l^)(2/,)3(^-:^). 



;r2r(3-^)|£(,.r^(l + i) 



X 1 + 



12c2 



exp<- 



2£„ 



3|£fe)l 



1 - 



12c2 



(75) 



In Eq. (j74ji an expansion on the Ipjc^ parameter has been used 
to clearly indicate the /p/c^-scaling. Without /^/c^-expansion 
the exponential factor of the differential ionization rate reads 



dw 
d^ 



exp ' 



2V3HV3 



V3S dp 



(1 + E^)\E(7js)\ £(/?,) 



/8 



E(ri,) V 3 ^ ■ ©2 
where H = -y/r^^T/2(V¥^T^^^, T = 1 - /^/c^, and 



(76) 




(77) 



+ ^^(S2 + 1)^((Vh^+i)h2- Vs^+i). 



Using the expression Eq. ( [76| ), the integration over momentum 
space can be done analytically. It yields for the total ionization 
rate: 



w - 



1 + 



161/p 
72c2 



exp 



4/. 



r(3 -^)^ ^ 

3£o\ 12c2j 



exp 



(78) 



In Fig. |4]the total ionization rate of Eq. (78 1 is compared with 
the ITM result of Ref.fS^]. The Coulomb-corrected relativistic 
SEA and the relativistic PPT yields results for the total ioniza- 
tion rate which are close, though they are not identical. In the 
next section we modify the Coulomb-corrected SEA with the 
aim to obtain ionization probabilities closer to the relativistic 
PPT result. 



III. MODIFIED COULOMB-CORRECTED RELATIVISTIC 
STRONG-FIELD APPROXIMATION 

In the standard SEA in the Goppert-Mayer gauge, the transi- 
tion matrix element is given by Eq. Q where 



Hi„, = r ■£(! -k-o;) 



(79) 




FIG. 4: The total ionization rate vs ionization potential for Eo/Ea = 
1 /32: (blue, short-dashed) in the relativistic Coulomb-corrected SFA 
via Eq. ( |78| l, (red, long-dashed) in the relativistic Coulomb-corrected 
modified SFA via Eq. ( |98[ l, and (black, solid) relativistic PPT from 
Ref. 132]. 



and the wave function of the initial state <pf\r, t) describes the 
free evolving atomic system, fulfilling the equation 



(/5,-//o)<^f'(r,f) =0, 



(80) 



with Ho^H- Hi„, = Ho ^ca-p+pc^ + V^'Xr). However, 
in general, as it is pointed out in 1551 l56ll . the SFA can be 
developed employing different partitions of the total Hamil- 
tonian which will result in different physical approximations. 
The main deficiency of the Coulomb-corrected SEA based on 
the standard partition indicated above, is that the influence of 
the laser field on the atomic bound dynamics is completely 
neglected. In this section we propose a modified version of 
the SEA employing another partition of the Hamiltonian in the 
SFA formalism. The main motivation is to use such a parti- 
tion where the laser field has a contribution to the bound state 
evolution. In particular, we use the following partition: 

H = Ho + Hi„„ (81) 
//() = ca ■\p-^{r I c\ + pc^ + V^'\r) (82) 

Hi„, = r-E(?7). (83) 

In this partition the transition matrix element reads 



MJ.'^-i \ dtd^rt}/)(r,t)Hi„,4>y(rj), 



(84) 



where the bound state dynamics is determined by the following 
time-dependent Dirac equation: 

id,4,f = {cor ■ [p - k (r ■ E(?7)) lc\ + (ic^ + V^'\r)\ 4>f\ (85) 

In this specific partition the initial bound state wave function 
* is the eigenstate of the instantaneous energy operator 

& = ca(p + A(/7)) + /3c^ (86) 



J, with A{ri) - -k(r ■ E(ri))/c in the Goppert-Mayer 
gauge, see Eq. ([9]l. Note that in the nonrelativistic limit the 
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standard SFA in the length gauge corresponds to the partition 
in which the initial bound state is the eigenstate of the energy 
operator ||55|| . This point provides another argument in favor 
of the applied partition in the relativistic case. 

The term incorporating the laser electric field in Eq. ( 85 i 
describes the spin dynamics of the bound electron before ion- 
ization and has to be considered in the further calculation. We 
are concerned by the spin dynamics in the bound state and by 
the dynamical Zeeman splitting of the bound state energies due 
to the spin interaction with the laser magnetic field. Therefore, 



in solving Eq. ( 85 i we will restrict the bound-bound transi- 



tions only to the transitions between dififerent spin-states at 
fixed quantum numbers {n, j, nij). The dipole approximation 
for the laser field will be adopted for the description of the 
bound state evolution: E(?7) « E(f), since the typical length- 
scale for the bound dynamics is much smaller than the laser 
wave length: rhIA ~ oj/kc ~ -/(Eo/EaXK/c) <k 1. When 
the initial spin-polarization is along the laser magnetic field 
l^^'^)) - \(p^g_J, then no spin transitions occur because the inter- 
action term Viat(t) = -(a ■ k)(r ■ E(f)) does not cause spin-flip 

(c) (c) 

{4>y:f\Vint{t)\4>g - 0. Accordingly, in this case we are looking 



for the solution of Eq. ( 85 1 in the form: 



where (p^^Jj,t) - 0^'^(r)e"'^±' is the ground state before 
switching on the laser field, i.e., an eigenstate of the atomic 
unperturbed Hamiltonian 

[ca-9+ pc" + y(^'(r)] = sJ^J.r). (88) 

The spin quantization axis for ip^^jx, t) states is along the laser 
magnetic field direction. Taking into account that 



(87) 



f 

*J — c 



./.'<<ilv..(0lO-^(i-il)' 



(89) 



where A(ri) = - J^^ E(ri')dT]', the wave function of the bound 
state will read 



_(f)exp 



2c \ 



2^ 
3c2 
21, 



1 - 



p_ 

3c2 



(90) 



The Zeeman splitting term in the phase of the wave function 
is of order of A/c ~ E()TcIc ~ k/c <& 1. It is small which 
arises from the fact that the perturbation term yint(f) couples 
the large/small spinor part of the initial state bispinor with the 
small/large spinor parts of the final state bispinor 

The explicit condition justifying the neglect of transitions to 
excited states can be given as follows. The transition prob- 
ability between the states « — > n' with energy diflTerence 
(^nn' ~ Ip (n and n' are the principal quantum numbers) is 
P„„i ~ Eorb/oL>„„' ~ Eorh/Ip, while the probability of a spin- 
transition in the state n is Pss' ~ E^rhTK {s and s' are the spin 
quantum numbers). Therefore, 



p , 

-* n IT 
P , 



1 



Eo 

E' 



(91) 



and the transitions to excited states are suppressed by a factor 

of Eo/Ea- 



With the wave functions of the bound state Eq. (90i, the 
differential and the total ionization rates can be calculated in 
the same way as in the previous section before. The structure 
of the ionization matrix element remains the same 



= -/m/,('?.)exp<j- 



[2Ms-c^+Ip)-pl] 
3|E(77o)|A 



3/2- 



(92) 



where 111/1(1],) is given by Eq. (59 1, in which the following re- 
placement should be done in order to account for the difference 
in the spinor operator of the interaction Hamiltonian in the 
modified SFA: 



( 1 - k ■ or) V/ 



1 + 



A(?7.5) - ad -H k ■ or) 
2^ 



(93) 



where v,- = v/ exp[+/A(?7)(l - 2Ip/3c^)] is the spinor part of the 
states given in Eq. ( [90| . The spinor operator in brackets in Eq. 
( 93 I comes from the spinor part of the Volkov wave function. 



Note that the phase in the exponential in Eq. ( 90 1 varies slowly 



compared with the contracted action S . Therefore, it does 
not modify the saddle point integration, but alters the pre- 
exponential term. 



y2ce^~ (2c + ips) 



^8c4 + 6c^pI 



(94) 



■Pi 



0^e-'^{2c + ipE) (8c5 - AQic^pE + Uc^p^^ 

12V2(8c4 + 6cV|+P^f' 
I3^e-'^{2c + ipE) {-2Sic^pl + 3cpl - 4-ipl) 



12V2(8c4 + 6c2p2 +p4^ 



3/2 



V2ce (2c - ipg) 



(95) 



ISc*+6c^pI+pI 
fe^ (8c5 + 40/cVe + 14c VI + 28ic^pl + 3cpl + 4i>|) 
12 V2(2c + ipE) (2c2 + pI) ^8c4 + 6cV| + pt 

and m+_ - m_+ - O. With this the differential ionization rate 
in the modified SFA yields: 



dw ^ ' -t- I2c2 24c« 



^)(2I,9i'-'^)co..p{'-^) 



.^r(3-?i^)|£(,.)|-^(i.i) 

41/„\ \-2{2A(s-c^+Ip)-plf^ 



X 1 + 



12c2 



exp 



3|E(7?,)|A 



(96) 
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and the total ionization rate is 



Tin 



exp 



41, 



(97) 



(21 pY 



■ exp 



12c2 



In Fig. |4j the total ionization probability calculated via the 
modified Coulomb-corrected relativistic SFA is compared with 
the result of the standard SFA, as well as with that of the rela- 
tivistic PPT. The conclusion can be drawn that both relativistic 
Coulomb-corrected SFA give slightly different results com- 
pared to the PPT for large Ip. The modified SFA is closer 
to the PPT result than the standard SFA. The total ionization 
probabilities in the relativistic Coulomb-corrected standard 
SFA is larger than the relativistic PPT by a factor smaller than 
1 + 3Ip/c^. 

The standard SFA yields larger ionization probabilities than 
the modified SFA which can be explained by a simple tunneling 
picture. In the first case the magnetic field acts on the electron 
only when it enters the barrier. In this case during the tunneling 
the electron kinetic energy will change on the amount of the 
interaction energy with the magnetic field /iZ?o = -sEqIIc {s - 
+ 1 for spin along the magnetic field or opposite), giving rise to 
a kinetic energy splitting for different spin orientations of the 
electron and, consequently, to different tunneling probabilities. 
Whereas, in the second case the magnetic field is also acting 
before tunneling and no kinetic energy splitting occurs. The 
ratios of the Keldysh-exponent for these two scenarios can then 
be expanded in the spin energy: 



exp 



2(2(7, +£o/2c))^' 
3£o 



-I- exp 



2(2(/p^£o/2c))'- 
3£o 



- modified 



2 exp 



2(2(7,,))'' 



3£|, 



1 + Iplc" 



(98) 



To decide which SFA partition is best suited to model the 
ionization process, a comparison with numerical simulations 
or experimental data is necessary. 

Further, we note that the Coulomb-correction term J dria ■ 
Vy'^^Vc that is neglected in both strong field approximations 
gives a correction factor of only 1 + (Ip/c^) V-EoZ-Ea, which is 
of the order of a few percents for the most extreme parameters 
(£o/£a = 1/16,Z = 90). 



into the relativistic regime. The applied Coulomb-corrected 
strong-field approximation incorporates the eikonal-Volkov 
wave function for the description of the electron continuum 
dynamics. The latter is derived in the WKB approximation 
taking into account the Coulomb field of the atomic core per- 
turbatively in the phase of the WKB wave function. In physical 
terms, the disturbance of the electron energy by the Coulomb 
field is assumed to be smaller with respect to the electron en- 
ergy in the laser field. The eikonal-Volkov wave function is 
applicable when the laser field is smaller than the atomic field 
Eq <k Ea and Ip < c^. 

We have derived an analytical expression for the ionization 
amplitude in a linearly polarized laser field within the Coulomb- 
corrected relativistic SFA when additionally smallness of the 
parameters of Ip/c^ <k and y «: 1 is used. A simple expression 
for the amplitude is obtained when using the Goppert-Mayer 
gauge. Moreover, in this gauge a Coulomb correction factor 
(ratio of the Coulomb corrected amplitude to the standard SFA 
one) coincides with that derived within the PPT theory. The 
differential and total ionization rates are calculated analytically. 
The calculated total ionization rate is slightly larger than the 
PPT rate at large ionization potentials. To improve the pre- 
dictions of the Coulomb-corrected relativistic SFA, we have 
proposed a modified SFA approach which is based on another 
partition of the total Hamiltonian. In this approach the SFA 
matrix element contains the eigenstate of the energy operator 
in the laser field as the wave function of the initial bound state. 
The modified SFA takes into account the dynamical Zeeman 
splitting of the bound state energy due to the spin interaction 
with the laser magnetic field (when the electron initial polar- 
ization is along the laser magnetic field) and the precession of 
the electron spin in the bound state (when the electron initial 
polarization is along the laser propagation direction). 

Our results show that the SFA technique allows the analyt- 
ical calculation of quantitatively correct differential and total 
ionization rates in the relativistic regime which takes into ac- 
count the impact of the Coulomb field of the atomic core as 
well as the electron spin dynamics in the bound state. This 
method can be viewed as an alternative to the PPT. In the fol- 
lowing paper of the sequel the Coulomb-corrected relativistic 
SFA will be used to investigate spin-resolved ionization prob- 
abilities. While for the total ionization rate the prediction of 
the standard SFA is close to that of the modified SFA, for spin 
effects their predictions are quite different. The next paper of 
the sequel will be devoted to this issue. 



IV. CONCLUSION 
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